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C*") . At hadron colliders, one of the most important channels for sparticle production is expected 

to be the gluino pair production. In the scenario where the sbottom is lighter than the gluino, 
gluinos can decay into sbottom and bottom quarks. Sbottoms can subsequently decay into 
bottom quarks and neutralinos. Hence, one expects a rich signature consisting of four b-jets 
and the missing transverse energy from neutralinos. To compute the cross section of this 
reaction, it is important to find tools and techniques suitable for the calculation of diagrams 



reactions. 



• with many particles in the final state. We present here a way to deal with such complicated 
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"^h ', 1 Introduction 

In the Minimal Supersymmetric extension of the Standard Model (MSSM), every funda- 
mental particle must have a supersymmetric partner with spin differing by 1/2 unit. When 
■ considering the conservation of R-parity, the lightest supersymmetric particle (LSP) must be 

stable, and supersymmetric particles are produced only in even number. If the LSP is also 
neutral and colorless, its signature in detectors will only be a large missing transverse energy. 



One often assumes the lightest neutralino to be the LSP. 



The main aim of present and future hadron colliders is the search for beyond Standard Model 
physics, and in particular Supersymmetry. Because of the high center of mass energies (TeV 
scale), the final states will be very complex, and multi-particle final state reactions will be pro- 
duced at a huge rate. Moreover, in order to have a good description of the physics at these 
scales, it is necessary to have precise theoretical predictions. 

In this work, we consider the production and decay of gluino pairs at hadron colliders, which 
appears as one of the most important channels for sparticle production. We investigate ways to 
achieve high precision calculations for multi-particle final state amplitudes in the MSSM. In par- 
ticular, while considering the complete gluino pair production and decay chain, one has a eight 
body amplitude which requires very complicated calculations. This amplitude is an interesting 
signature at Tevatron. Furthermore, it is a good example to develop a computation technology 
able to treat complex reactions. 

In section 2 we give an overview of the gluino pair production. Section 3 describes the gluino 
decay modes and presents the topologies. The way we performed the calculation of the corre- 
sponding amplitudes is presented in section 4. Finally, conclusion and perspectives are given in 
section 5. 
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2 Gluino pair production 



Gluinos and squarks are presently actively sear ched at Tevatron. Lower bounds have already 
been set up for the masses of gluinos and squarks ^^ l In particular, if gluinos and squarks have 
the same mass, the lower mass limit is 310 GeV. In the opposite case, if the gluino (squark) is 
heavier, the lower squark (gluino) mass limit is then 240 GeV. 

The main modes for gluino pair production at hadron colliders are, in leading order of the per- 
turbative expansion: 
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The relative weights of each reaction is strongly correlated to the masses of squarks and gluinos^. 
Due to the presence of an extra intermediate sparticle in the t-channel, for both cases the s- 
channel is privileged. 

3 Gluino decay 

In the MSSM, for the third generation of quarks, a strong mixing can appear for the superpartner 
masses, which depends on the parameters of the theory. In particular, the masses of the two 
sbottoms read: 
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where the dependence on tan (3 (ratio of the VEV's of the two Higgs fields) and Aj, (sbottom- 
Higgs trilinear coupling) appears explicitly. Thus, for a large value of tan /?, one observes a 
high mass splitting for the two sbottoms. This can lead to a rather low mass for the lighter 
sbottom, so that the g — > b + b decay would be possible. The sbottom can then decay through 
b^b + xl 

The complete reaction to be considered is therefore qq / g g — > b x± b b Xi b. We select the 
topologies represented schematically in figure Q which can be experimentally identified through 
kinematic cuts. The experimental signature will then be the missing transverse energy of two 
undetected neutralinos and the presence of four b-jets. To compute this amplitude, one has to 
consider all the leg permutations. This leads to 16 diagrams for each of the s-channels (a, c) 
and 32 diagrams for the t-channels (b, d), so that 48 diagrams have to be computed for each 
partonic subprocess. 



4 Cross section calculations 



In usual methods for computing the amplitude of diagrams with many particles in the final state, 
one generally splits the processes into several different subprocesses, and obtains the total cross 



(a) (b) (c) (d) 

Figure 1: Schematic representation of the topologies corresponding to qq / g g — > b \i b b xi b. (a) and (b) 
correspond to the gluino pair production from qq annihilation, for respectively s- and t-channels. (c) and (d) 
correspond to the gluino pair production from gg, for respectively s- and t-channels. 



section while considering the branching ratio of each subprocess. The computations are then 
relatively easy, but one loses the information about particle spins and one obtains approximate 
results. 

On the contrary, the direct calculation is an exact method, which consists in calculating the 
total amplitude of a process and integrating it numerically over the phase space associated to 
the final state particles. This method is unfortunately difficult to follow, especially when the 
number of diagrams corresponding to the studied process is large or/and if the number of par- 
ticles in the final state is high. To achieve those computations, sev eral automatic codes (such 
as CompHEpEI GRACE ^, MadGraph^, AlpgenEH, FormCalcUl, ...) exist, but only a few 
are able to deal with multiparticle processes in the MSSM. In particular, CompHEP is almost 
able to treat completely 2 — ► 6 MSSM amplitudes. Unfortunately, in this case, the number of 
generated terms in the output files is so high that the results cannot be used, neither analytically 
nor numerically. 

We computed this amplitude with Formal and Mathematical^, using a modified version of the 
FormCalc package. 

In its standard version, FormCalc can deal with 2 — > 3 amplitudes within the MSSM. Prac- 
tically, using first the FeynArtsEl part of the package one can generate Feynman diagrams 
and write the corresponding Feynman amplitudes. Up to this step, it is possible to generate 
amplitudes with as many particles in the final state as wanted, but subroutines to go further 
do not exist for more than three particles. In the next (internal) steps, FormCalc uses Form to 
perform analytical calculations of the amplitudes, and then generates a Fortran code containing 
the squared matrix elements. Finally, it computes the cross-sections of 1 — ► 2, 2 — > 2, and 
2 — ► 3 processes. 

In order to compute our 2 — ► 6 amplitude, we had to modify the internal files of Form- 
Calc which are written in Form. Those modifications are suitable for any 2 — ► 6 amplitude, 
and can be easily extended to even more complicated cases. The generated Fortran code for 
the squared matrix element s is then relatively concise and fast thanks to the use of the Weyl- 
van-der-Waerden formalism^! Instead of elaborating new kinematics modules for FormCalc 
enabling to calculate numerically the cross section, we chose rather to insert our results into 

pythiaED. 

Complete results will be presented in a future publication. 



5 Conclusion and perspectives 



The gluino pair production is known to be an important channel for sparticle production at 
hadron colliders. In this note, we considered the gluino pair production from qq and gg initial 
states, and their decay into sbottom and bottom quarks in the MSSM. The whole reaction is 
then qq I gg — ► b Xi b b Xi b. The direct calculation of this amplitude has been performed 
while modifying the FormCalc package for Mathematica. This study allowed us to get adequate 
tools and techniques for the direct calculation of 2 — > 6 amplitudes. 

The next step, in collaboration with the DO group at IPNL, consists in a comparison of the 
theoretical results and the experimental data, in order to derive constraints on the gluino and 
sbottom masses. Moreover, it will be interesting to compare the results of the direct calculation 
to those of the usual method used in PYTHIA (branching ratio method) and to check whether 
this simplest method is reliable for complicated cases with many particles in the final state and 
asymmetric events. 

The technology employed for these calculations can hopefully be used for other amplitudes 
or for even more complicated cases. Thus, it will be very useful and interesting to use these 
tools to make precise predictions for Tevatron and LHC. 
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